Impurity mediated nucleation in hexadecane-in-water emulsions 
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We report detailed nucleation studies on the liquid-to-solid transition of hexadecane using nearly 
monodisperse hexadecane-in-water emulsions. A careful consideration of the kinetics of isothermal 
and nonisothermal freezing show deviations from predictions of classical nucleation theory, if one 
assumes that the emulsion droplet population is homogeneous. Similar deviations have been observed 
previously || . As an explanation, we propose a novel argument based on the dynamic generation of 
droplet heterogeneity mediated by mobile impurities. This proposal is in excellent agreement with 
existing data. 



PACS numbers: 82.60.Nh, 64.60. Qb,44.60.+k,82. 70. Kj 



I. INTRODUCTION 

The use of micron-sized emulsified droplets for nucle- 
ation studies of the liquid-to-crystal transition is a well 
established technique. Vonnegut Q pioneered this con- 
cept in the first half of this century and others have used 
it to study a variety of materials The advantage 

of using a large ensemble of independent nucleation sites 
to measure the stochastic process of nucleation is obvi- 
ous. An additional advantage is that the effects of crys- 
tal growth are virtually eliminated. On the time scale of 
the nucleation measurements, the growth time for each 
droplet is instantaneous. Furthermore, because exper- 
iments are conducted in a range where the probability 
of two nucleation events per drop is small, it is a simple 
matter to equate the total number of crystallized droplets 
with the number of nucleation events. 

The role of impurities on nucleation in emulsion stud- 
ies is not obvious. Clearly, one possible role is to act 
as heterogeneous nucleation sites. Both Turnbull Q and 
Perepezko || have discussed this. A signature of het- 
erogeneous nucleation is a small difference between the 
melting temperature and the onset of nucleation. For 
example, in mercury, Turnbull showed that changes in 
surfactants could increase the undercooling from 5 °C to 
60 °C, with the smaller values attributed to the effects of 
heterogeneous nucleants either on the surface or within 
the volume of the droplets. It is well known that emul- 
sification tends to increase the undercooling over that of 
bulk liquids, but a simple calculation shows that for this 
to result from isolation of the heterogeneous catalysts 
they must be present at extremely low levels (cf. Sec. 
0). This seems to suggest that most impurities are not 
effective nucleating agents ||. 

A second role of impurities is to lower the melting 
point. Perepezko || has shown in his studies of metals 
that as the liquidus temperature drops due to alloying, 
there is a corresponding decrease in the nucleation tem- 
perature. Hence an approximately constant value of the 
undercooling is observed. These effects can be large, for 
example in the lead-antimony system the liquidus and 
nucleation temperatures both decline by about 75 °C as 
the antimony content increases up to 17.7 atomic per- 



cent. The effect of a typical impurity level (mole fraction 
ph 0.01) will of course be much smaller. However, the 
nucleation rate is strongly temperature dependent. For 
example, in n-alkanes (C n H2 n +2, henceforth abbreviated 
On) , the nucleation rate changes by a factor of 5000 per 
°C (a). We show in this work that impurity levels of even 
a few percent cause a significant change in nucleation be- 
havior. This influence must be accounted for to extract 
an accurate value of the nucleation rate. 

Size of the emulsion droplets also plays a key role in 
nucleation studies. In homogeneous nucleation, the nu- 
cleation rate is proportional to the volume of the droplets. 
Typically, the determination of the size distribution for 
the emulsions is a large source of error in nucleation rate 
measurements ||. Advances in emulsion synthesis tech- 
niques now makes it possible to create nearly uniformly 
sized emulsion particles. The average deviation in diam- 
eter from the mean size is only from 10 to 15% [0-^]. By 
using such emulsions with narrow size distributions, the 
determination of the rate constant for any given emul- 
sion radius is considerably improved. Also, this makes 
it possible to better test the predicted volume scaling of 
the nucleation rate. 

In this paper, we report the results of a refined exper- 
imental and theoretical investigation of the nucleation 
rate in emulsified hexadecane (C16). The quality of ex- 
perimental data has been significantly improved by using 
nearly monodisperse emulsions in a well controlled ther- 
mal environment and by using x-ray scattering to accu- 
rately monitor the volume of nucleated droplets during 
crystallization. Data were obtained for both fixed and 
linearly increasing undercooling as a function of time, 
henceforth referred to as isothermal nucleation and lin- 
ear cooling, respectively. Thermodynamic melting curves 
have also been obtained for the samples in order to as- 
sess the influence of impurities introduced during the 
emulsification procedure. The theoretical analysis has 
concurrently been refined to account for the entire time- 
dependence of the solid fraction as a function of time, 
rather than just matching some characteristic time. This 
careful analysis revealed deviations from a simple sce- 
nario involving a uniform ensemble of independent nu- 
cleation events. By systematically relaxing the assump- 
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tions in this description, we have concluded that the most 
probable cause of the observed deviations is a novel mech- 
anism which involves the transport of impurities expelled 
from nucleated droplets to the remaining liquid droplets. 
This transport subsequently increases the impurity con- 
centration of these remaining droplets and reduces their 
nucleation rate by lowering their thermodynamic melting 
point T m . This mechanism can completely account for 
the observed behavior and can be further explored with 
additional experiments. 

The remainder of this paper is organized as follows: In 
Sec. ||_ we review previous studies on alkane nucleation. 
Sec. Ill provides details of t he experimental techniques 
and procedures used. Sec. IV reviews classical nucle- 



ation theory and its predictions based on the assump- 
tion of uniform nucleation, i.e., independent nucleation 
events in a homogeneous ensemble of droplets. Finding 
that this theory is unable to account for the observed 
behavior of the system, we consider possible extensions 
of this theory in Sec. [V| Here we explore two scenar- 
ios for the generation of droplet heterogeneity, which is 
necessary to explain the data. These include the possibil- 
ity of some hxed-in-time heterogeneity generated by the 
emulsification process, such as a distribution of droplet 
sizes or impurity concentrations, and the dynamic gener- 
ation heterogeneity in the droplet population due to the 
nucleation process itself. Based on the reanalysis of the 
data, we present our conclusions in Sec. VI and suggest 
a number of experiments to further test the validity of 
the proposed mechanism, as well as avenues for future 
theoretical consideration. 



II. PREVIOUS STUDIES OF ALKANE 
NUCLEATION 

Four groups have studied alkane nucleation through 
the use of emulsion samples. The earliest work is from 
Turnbull and Cormia § who studied C16, C17, C18, C24 
and C32. Theirs is the first evidence that alkane nucle- 
ation is unusual. First, they noted that there seemed 
to be an unusual spread in the melting temperatures. 
This was characterized as "sharp" and "broad" melting 
fractions. The amount of each varied from sample to 
sample, even for the same chain length. Because they 
performed isothermal nucleation studies they were sensi- 
tive to this spread in melting temperature. To analyze 
the nucleation behavior in this situation they focused on 
the early-time data with transformed fraction n < 0.5. 
This narrowed focus was meant to isolate the behavior of 
the sharp-melting, majority phase. 

The second anomaly in alkane nucleation is the ease 
with which the alkanes nucleate. Stated in terms of re- 
duced undercooling, AT r = (T m — T/v)/T m , where Tjq 
is the point where the nucleation rate becomes signifi- 
cant and T m is the thermodynamic melting temperature, 
AT r for the alkanes is about 0.05 whereas that for other 



materials is from 0.2 to 0.5 Turnbull and Cormia's 
analysis of the nucleation behavior in terms of the classi- 
cal nucleation model showed that the nucleation barrier 
is small (9.64 mJ/m 2 for C18; detailed kinetics measure- 
ments were not performed for C16) but that the pre- 
exponential factor is in rough agreement with that cal- 
culated from classical nucleation theory (experimental 



value 



10' 



37.35±2 



m 3 s for C18). The small barrier 



accounts for the small undercooling temperature. The 
agreement of the pre-exponential with that from classi- 
cal nucleation theory is itself unusual. Other materials 
typically exhibit values several orders of magnitude larger 
in value than the classical value [Q . 

Two subsequent studies on alkane nucleation explored 
the effect of changing carbon number 0JT|. In both 
studies continuous cooling data were used to calculate 
the crystal-liquid surface tension and the pre-exponential 
factor. Both studies found that as the carbon number is 
reduced below about np = 15, there is an increase in the 
barrier. The pre-exponential values were not precisely 
determined by these studies. Re-analysis in terms of the 
Negentropic Model |l0| led to a value of 9 mJ/m 2 for the 
barrier of C16, and to the suggestion that the anomalous 
barrier height disappears as fip — > 1 and as nc — ► oo, 
and that its origin is in partial alignment of the alkane 
chains. 

The fourth group [T^-[l^] to study nucleation in emul- 
sions focused on the behavior of C16. They utilized ul- 
trasound transmission to measure the proportion of liq- 
uid to solid in an emulsion sample. Both stepwise cooling 
and isothermal hold experiments were used. Their results 
exhibit the typical ~ 14 — 15°C undercooling found by 
other workers, showing that their samples behave as ex- 
pected during stepwise cooling. However, the real focus 
of these studies was the interaction of liquid and crys- 
talline droplets during isothermal hold experiments p3[ . 
At an undercooling of ~ 10 °C, where the nucleation rate 
for the C16 liquid droplets alone is quite small, equal 
volume mixtures of liquid and solid C16 droplets were 
studied. The solid C16 drops were created by an initial 
deep undercooling and then physically mixed together 
with the as-yet unfrozen sample. The authors found that 
the solid C16 particles accelerate nucleation. The parti- 
cle size distribution before and after the experiment was 
similar; therefore they ruled out Ostwald ripening as the 
source of the acceleration. No attempt was made to fit 
this data to a rate law. 

In a set of related experiments the type and concen- 
tration of the surfactant was varied |l4|]. The presence 
of crystalline C16 droplets accelerated the nucleation by 
varying amounts depending on the surfactant type. It 
was also found that the rate of nucleation for the mixed 
solid/liquid emulsion increased as the amount of surfac- 
tant increased. The authors suggest that inter-particle 
collisions are responsible for the accelerated nucleation. 

As we discuss in Sec. IV, we observe an effect opposite 
to that seen by these authors. As our sample crystallizes, 
the remaining crystallization is more difficult. However, 
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FIG. 1. Left: A schematic of the x-ray cell loaded with 
an emulsion sample. Right: The powder pattern of the C16 
emulsion sample in the triclinic (solid) phase. Numbers in- 
dicate the Miller indices for each peak. The arrow denotes 
the detector position, which monitors the (013) peak to de- 
termine the amount of solid in the sample during the kinetics 
measurements. 



Melting Temperature: T m o 


291.32 K 


Entropy of Fusion: AS" 


6.28 x lO^Jm^KT 1 


Mass Density: p 


773.4kg m -3 


Molar Weight: M w 


0.22643kg moh 1 


Viscosity at T m o: i) 


3.484 x 10 _;i Pa s 



TABLE I. Material Properties of Hexadecane |[18 



Dickinson et al. |l6| ] also report an experiment where this 
effect is seen. In this work they hold their emulsion at 
a undercooling of about 15 °C and follow the percent 
transformed with time. Like us, they observe an initially 
rapid nucleation which then slows. They do not attempt 
to fit this to a rate law, but it is clear that a single ex- 
ponential cannot describe their data. Also similar to our 
findings they report that the melting behavior extends 
over a large temperature regime. They report melting 
ranges of about 2 °C for all their samples, and in one 
instance [ ft6| they report the onset of melting more than 
10 °C b elow the bulk melting temperature. As we show 
in Sec. V A, this is a signature of the impurity effects 
which can affect nucleation. 



III. EXPERIMENTAL PROCEDURE 

Nucleation of the liquid to crystal transition was mea- 
sured in droplets of C16 of an oil-in-water emulsion. The 
observed crystal phase is the thermodynamically stable 
triclinic phase. A schematic of the experimental setup is 
shown in Fig. [l]. A reference table of relevant material 
properties for C16 is given in Table |. 

One fim diameter emulsions of C16 in water with the 



surfactant sodium dodecyl sulfate (SDS) were prepared 
using the fractionation method of Bibette [Jl7| . The dis- 
tribution of droplet sizes was less than ±10% in diameter 
as determined by light scattering and optical microscopy. 

Below its Krafft point at 283 K Q, SDS precipi- 
tates out of solution as a crystalline solid in equilib- 
rium with a small concentration of dissolved monomers. 
In order to avoid this effect, a co-surfactant was added 
to supress the Krafft point j2(J. The chosen sur- 
factant, a sulfated polyoxyethylenated alkyl alcohol 
(refered to here as EO), has the chemical formula 
CH 3 (CH 2 )i 1 _i5(OCH2CH 2 ) 9 OS0 3 Na. An aqueous so- 
lution of EO was added to the initial emulsion sample to 
achieve approximately equal concentrations of SDS and 
EO in the range of 0.5 - 1 wt%. This new solution was 
allowed to equilibrate for more than 24 hours before use, 
during which time the emulsion droplets had flocculated 
together to form a white cream floating on the aqueous 
phase. The cream of the emulsion sample with both SDS 
and EO co-surfactants was loaded directly into the x-ray 
cell. 

The emulsion sample was loaded into a Viton O-ring 
sandwiched between two Be windows, held between two 
copper plates (see Fig Yu). The sample temperature was 
read with thermistors (±0.2 K accuracy) placed in the 
copper less than 3.5 mm from the sample. The sam- 
ple was cooled on an Air Products Displex cold finger 
with a closed-cycle He refrigeration system. Heating wire 
wrapped around the base of the cold finger provided tem- 
perature control. As the cold finger was operated in vac- 
uum, an air-tight copper shield with kapton windows was 
placed directly around the copper sample mount inside 
the vacuum shroud of the Displex to surround the sam- 
ple with gas and therefore reduce temperature gradients. 
The temperature was controlled with a LakeShore 340 
Temperature Controller. The control temperature sen- 
sor (Si diode) was embedded in the cold finger directly 
above the copper sample mount. Using this setup, it was 
possible to cool the sample from room temperature down 
to 277 K in less than 250 s without any overshoot of the 
temperature (Fig. ||). This careful attention to temper- 
ature gradients overcame problems encountered in the 
initial studies |plfl . 

The nucleation rate was determined by measuring the 
change in intensity of the x-ray diffraction as a function 
of time. X-ray scattering was performed using a Rigaku 
18 kW rotating anode generator and an x-ray wavelength 
of 1.542 A. Scattered x-rays were detected with a Bi- 
cron detector on a Huber 4-circle spectrometer. Scat- 
tering from the solidified 1 fim dia. droplets resulted in 
a powder diffraction pattern of the triclinic phase. The 
(013) Bragg peak of the triclinic phase was selected to 
be monitored (Fig [l]) due to the intensity of the peak 
and the low scattering from the liquid phase at that an- 
gle. For the isothermal experiments, the detector slits 
were set very wide to increase the signal intensity and to 
encompass slight differences in the peak position due to 
thermal expansion when measuring at different tempera- 
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FIG. 2. The temperature of the sample and measured x-ray 
intensity as a function of time over the course of a typical 
isothermal experiment. The sample was first held above T m to 
ensure that all droplets were liquid (1). The sample was then 
quenched to (2) and held at (3) a fixed temperature T, and 
the rise in scattering intensity was measured. A subsequent 
quench to a lower temperature (4) was used to determine the 
scattering from a fully crystallized sample. The sample was 
then heated back above T m (5). 

tures. During the constant cooling rate and the melting 
experiments, the slits were narrowed and the detector 
was continually shifted to track the peak during thermal 
expansion. 

Three types of experiments were conducted: crystal- 
lization at a constant undercooling (isothermal nucle- 
ation), crystallization while the temperature was lowered 
linearly in time (linear cooling), and melting at a slow 
heating rate. For the isothermal experiments, the sam- 
ple was rapidly quenched from room temperature down 
to the desired undercooling and held at that temperature 
while the increase in x-ray scattering was followed as a 
function of time (see Fig. |^). In order to determine the 
total scattering signal when the entire droplet popula- 
tion was frozen, each experiment was followed by a deep 
quench to a lower temperature (273.8 K) that froze all 
droplets. All reported isothermal data is normalized to 
the final intensity value of the subsequent deep quench. 
In addition, all quenches were repeated multiple times 
with identical results. 

In the linear cooling experiments, the rate of temper- 
ature change was between —0.01 and —0.002 K/min. 
Faster rates were avoided because the latent heat re- 
leased by the crystallizing droplets could not be dissi- 
pated quickly enough and resulted in a smearing of the 
transformation over a larger temperature range. 

Melting was done in steps of 0.25 to 1 K. The sam- 
ple was heated slowly between steps and was allowed to 
equilibrate for approximately 1 hour at each step. This 
equilibration was observed by the leveling out of the in- 
tensity subsequent to each temperature step. 



IV. UNIFORM NUCLEATION 

In this section, we interpret the results of the experi- 
ments using classical nucleation theory, with the assump- 
tions that (i) all droplets in the emulsions are identical, 
and that (ii) nucleation events are statistically indepen- 
dent in each droplet. We refer to this set of assumptions 
as uniform nucleation. 

In classical nucleation theory (Jj2^], the formation of 
the thermodynamically stable crystalline phase in an un- 
dercooled, metastable liquid phase is controlled by the 
local free energy barrier towards solidification. The free 
energy cost, AG(R), associated with the creation of a 
solid sphere of radius R is given by 

A<jr 

AG = -—R 3 Ag + 47ri? 2 7, (1) 
o 

where Ag is the free energy difference per volume be- 
tween the metastable liquid and the stable solid phase, 
and 7 is the surface tension between the two phases. 
For small solid radii, R, the second term dominates the 
first and the formation of the more stable solid phase 
requires an increase in free energy. The nucleation rate 
is controlled by the free energy barrier associated with 
the formation of a solid sphere at a critical radius of 
R c = 2j/Ag, since further solidification reduces the 
free energy of the system. For small undercooling, the 
thermodynamic driving force Ag = AS AT, where AS 
is the entropy of fusion of the phase transition and 
AT = T m — T is the degree of undercooling of the sam- 
ple at temperature T [J23| . We neglect the dynamics of 
the growth and simply assume that upon formation of 
a solid phase of critical radius, the entire droplet freezes 
instantaneously. 

The probability of spontaneously creating a solid with 
the critical radius is simply found from Eq.(|l|), leading 
to a nucleation rate 

This rate is proportional to the volume, V, of available 
undercooled fluid, the Boltzmann factor associated with 
the free energy barrier for creating a critical nucleus as 
determined from Eq.(|l]), and an attempt frequency per 
unit volume, A, for creating such a critical nucleus, which 
is given by [Q 

Here, C, lu, a, d and D = (k B T) / (3irr]d) are, respec- 
tively, a numeric constant C = 1.65, the monomeric vol- 
ume, molecular number density, molecular diameter and 
the Stokes-Einstein diffusion constant of C16 in its liquid 
phase. For notational simplicity we define 



so that the nucleation rate given by Eq. (0) becomes 

SI 

v Q = AV exp 



(AT) 1 



(5) 



A. Isothermal Nucleation 

At a fixed level of undercooling, the rate at which emul- 
sion droplets freeze is time independent. Let n(t) repre- 
sent the fraction of solid droplets at time t, following a 
rapid quench to a given undercooling AT at t = 0. Then 
n(t) obeys the differential equation 



dn 
~dt 



(1 - n) v Q , 



(6) 



with initial condition: n(0) = 0. This differential equa- 
tion may be integrated to give 



n(i) 



1 



(7) 



so the fraction of liquid droplets decreases exponentially 
in time. 

It should be noted that the rate v a depends strongly 
on both the surface tension 7 of the liquid-solid interface 
and the degree of undercooling AT [cf. Eqs.(Q)-(]|)]. This 
strong dependence allows an extremely precise determi- 
nation of 7 provided that the data can be fit to a single 
exponential. On the other hand, Fig. |3|, showing a rep- 
resentative data set from an isothermal freezing experi- 
ment, demonstrates that such a simplistic interpretation 
of the data is inapplicable. The freezing of remaining 
liquid droplets at later times proceeds far slower than 
suggested by an extrapolation from the initial decrease 
of the fraction of liquid droplets, (1 — n), shown by the 
solid line in Fig. ||. 



B. Linear Cooling 

We now consider the case where the temperature, 
rather than being held constant, decreases linearly in 
time. The thermodynamic driving force thus increases 
in time and the exponential decay of the liquid fraction 
seen in the solution of Eq. (g) is changed to a relatively 
sudden drop at an undercooling determined by the cool- 
ing rate (see Fig. |^). 

Suppose that the temperature of the sample is de- 
creased at a rate given by A. If we ignore the weak 
temperature dependence in the pre-exponential factor in 
Eq. (Q), we find that the time- dependent nucleation rate 
is now given by 



v(t) = AV exp 
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FIG. 3. The fraction of liquid emulsion droplets as a func- 
tion of time at constant undercooling of 15.3 K as measured 
by the x-ray intensity at a Bragg peak of the crystalline C16 
phase. If the nucleation rate had been independent of time, 
this process should have yielded a simple exponential depen- 
dence on time. The solid line represents the expected behav- 
ior based on the nucleation rate at early times (t < 100 s.) 
The data deviates significantly from this simple model at later 
times. 



The differential equation describing the time evolution of 
n(t) reads 



dn 
~dt 



(1 - n) AVe~T*&, 



(9) 



where the initial time, t = 0, is chosen to occur when the 
sample is at T m . 

Following a change of independent variable to z = 
S7/ (Xt) 2 , the solution of Eq. (|ij) is given by 



n(t) = 1 — exp 



2X~ 



1 SI 



2' (xty 



(10) 



where we have introduced the incomplete Gamma func- 
tion pi 



1 

"2' Z 



e-'S-^dz. 



(11) 



For z ^> 1, this function goes to zero as z 3 / 2 exp(— 2 
so at short times n(t) behaves as 



n(t) 



AVX 2 
2S1 



t 3 e~^, t < X/VQ. 



(12) 



Thus, at t = the function has an essential singular- 
ity and the initial increase in the solid fraction n(i ) is 
extremely slow. For z <C 1, the integral in Eq. (|l l|) is 
dominated by the divergence of the integrand at z = 
and is approximately given by 2z~ x ^ 2 , so at long times 



(13) 



n(t) - 1 - exp(-AVt), t > A/VS1 
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FIG. 4. Solid fraction as a function of temperature at a 
constant cooling rate, A = 4.9x10 °C/min, plotted together 
with the theoretical curve given by Eq. (|l"o|). Note that time 
evolves from right to left. 



as expected since the free energy barrier disappears. 
However, since the attempt frequency AV is typically 
very large, the nucleation process is largely finished long 
before the free energy barrier becomes negligible. There- 
fore, this asymptotic behavior is not observable experi- 
mentally. The complete solution is plotted in Fig. |^ to- 
gether with the experimental data. 

The plot shown in Fig. | of Eq. @ was made using 
the pre-exponential facto r taken from a fit to the isother- 
mal data (see Sec. VB ) and using the surface tension 
7 as a fitting parameter. The pre-exponential factor A 
used in this fit was calculated at T m (A at the temper- 
ature of transformation is less than 3% smaller than A 
at T m ). We find that the best fit value of this surface 
tension is 7 = 9.9 ± 0.2 mJ/m 2 . It should be noted that, 
to a good approximation, shifting the value of the sur- 
face tension 7 simply translates the theoretical curve in 
time (i.e. temperature). Therefore, 7 can in principle be 
determined with great precision, since the time at which 
abrupt change in nit) occurs is exponentially sensitive to 
its value. The principal source of the uncertainty is the 
accuracy with which the magnitude of the undercooling 
was measured. This result for the surface tension is in 
agreement with the previous work of Turnbull et al. for 
C18 @ and other groups for C16 @@. 

Finally, we note that the agreement between the linear 
cooling data and the prediction of Eq. (|l(]) is not exact. 
In particular the slope of the data in the transition region 
(3.2 - 2.6 °C in Fig. |) is clearly smaller than that of the 
theoretical prediction. This behavior is consistent with 
the observed slowdown of the nucleation rate in isother- 
mal experiments. In order to address these discrepancies, 
we proceed in the next section to systematically relax the 
assumptions made in the uniform nucleation hypothesis. 



V. NONUNIFORM NUCLEATION 

Applying the results of uniform nucleation theory to 
both the isothermal and linear cooling experiments has 
led to substantial disagreements between theory and ex- 
periment in the former and more subtle deviations in the 
latter. The departure from simple exponential behavior 
in the isothermal quench experiments has been noted by 
previous workers in the field Faced with this ob- 

servation, one is forced to consider one of two general 
possibilities: 

(i) The emulsion droplets are not homogeneous, i.e., 
there is a "quenched-in" dispersion in some relevant prop- 
erty, such as a distribution of sizes or impurity concen- 
trations, leading to mute-exponential decay of n(t) in the 
isothermal experiments, or, 

(ii) Nucleation events are not statistically independent, 
i.e., the nucleation rate of a given emulsion droplet de- 
pends on the state of the other droplets, such as the 
current fraction of solid droplets in the emulsion, pos- 
sibly due to an impurity mediated interaction between 
the droplets. This sort of interaction might lead to 
a history dependent nucleation rate and result in non- 
exponential freezing curves. We re-emphasize that the 
non-exponential character of the curve (Fig. ||) does not 
reflect the finite growth rate of the solid as this rate 
is extremely fast compared to the data acquisition rate 
p6| , p7| . Similarily we may not attribute the deviation 
from a simple exponential to latent heat effects (rise in 
temperature due to rapid release of heat from crystal- 
lization) since the magnitude of the deviation does not 
diminish with smaller uncoolings and hence, lower nucle- 
ation rates. In addition, given that the isothermal curves 
can be reproduced repeatedly for the same loading of 
a given sample, the non-exponential character does not 
arise from coarsening of the size distribution of the emul- 
sion sample. The main difference between mechanisms 
(i) and (ii) is that in the former scenario, although there 
is a distribution of nucleation rates, it does not depend 
on the state of the droplets, whereas the latter scenario 
suggests that nucleation events influence the subsequent 
evolution of the system. 

Let us start by considering the first possibility, the ex- 
istence of some fixed in time distribution of nucleation 
rates. If there is a population of droplets with differ- 
ent volumes, Eq. ((7J) can easily be generalized to account 
for this effect by integrating contributions to the freezing 
curve of droplets of a given volume V over the volume dis- 
tribution V(V). The resulting isothermal freezing curve 
is then [cf. Eqs.(§), (g)] 



n{t) = 1- J dVV(V) exp (~AVe- n/(AT) '' 



(14) 



We take the droplet volume distribution to be log normal 
with characteristic volume Vq and variance a 1 : 28 



G 



V{V) 



V V2~t 



: exp 



2a 2 



(15) 



and find that the best fit to the experimental data is ob- 
tained for a « 1.0. The resulting best fit to the data 
is poor; also, the emulsion is actually known to have a 
much narrower volume distribution, based on the method 
of emulsion preparation and our own optical measure- 
ments p9[ . Furthermore, isothermal experiments at dif- 
ferent undercoolings cannot be accounted for by a single 
set of fitting parameters. Thus, we do not believe that 
the observed non-exponential saturation of the isother- 
mal freezing curves can be attributed to the polydisper- 
sity of the emulsion sample. 

Along a similar vein, one may postulate that there 
exists some fixed distribution of melting temperatures 
in the emulsion droplet population, possibly due to a 
distribution of impurity concentrations in the emulsion 
droplets. Here one must assume that the impurities in 
the droplets are highly insoluble in water, in order to 
prevent the equalization of impurity concentration via 
diffusion through the continuous phase. Problems very 
similar to the case of size distribution arise here as well. 
Although good fits to individual isothermal nucleation 
curves can be obtained using this modification, exper- 
imental runs at different undercoolings cannot be ac- 
counted for using the same fixed impurity concentra- 
tion distribution. Furthermore, in order to fit a given 
data set one is once again forced to assume a distribu- 
tion of impurity concentrations that is too broad (with 
standard deviation comparable to the mean). The av- 
erage number of impurities in a droplet of volume V is 
Ni = cpVNa/M w , where c is the mole fraction of im- 
purities, Na is Avogadro's number, and p, M w are the 
density and molar weight of C16, respectively. Substi- 
tuting material values from Table || and typical values 
c w 0.01 and V « 5 x 10" 19 m 3 , yields N « 10 7 . Thus, 
at this impurity level, the central limit theorem suggests 
a much narrower distribution. 

Nevertheless, the role of impurities in this system is 
not inconsequential. While it is true that we may dis- 
pense with a time-independent distribution of impurities 
in order to account for the experiments, a time-dependent 
change of the impurity concentration driven by the freez- 
ing process can account for the non-exponential behavior 
of the isothermal freezing data. Before we discuss this 
possibility in detail, we consider the effect of impurities 
on the melting temperature of the C16 droplets. 



A. Impurity Effects 

In a liquid, the presence of impurities that are rela- 
tively insoluble in its solid reduces its equilibrium freez- 
ing point |3^] and broadens its melting curve. Figure ^ 
shows the behavior for both the emulsion and the bulk 
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FIG. 5. Melting curves obtained for bulk and emulsified 
(d = 1.0 [im) C16 samples used in the nucleation study, 
shown as solid fraction, n, versus temperature. Solid lines 
correspond to best fits to Eq. ( |l8| ) with parameters quoted in 
the legend. 



C16 that was used in its preparation. Note that each 
emulsion droplet (or the entire bulk C16 sample) will 
have a partially molten fraction, described by x in the 
following discussion. 

For molar impurity fractions c«l, the impure liquid 
will be in equilibrium with its pure solid at a suppressed 
temperature [BOl 



T m (c) = T, 



mO 



(16) 



where T m o is the melting point of pure C16. A point on 
the melting curve at temperature T reflects coexistence 
between the liquid fraction x with impurity concentra- 
tion c/x and the pure solid, since the total number of 
impurities are fixed and the impurities are assumed to be 
insoluble in the solid. Thus, the temperature at which 
the liquid fraction is x is given by 



T(x) = T m (c/x) = T m0 - ac/x. 



(17) 



Here, a = k^T/ AS w 13.2 K, obtained from the material 
parameters of C16 in Table |. Inverting this equation 
gives the liquid fraction function of T : 



, T < T m0 - ac, 



x(c;T) = 




(18) 



A single-parameter fit to the bulk melting curve gives 
excellent agreement over the entire temperature range 
(see Fig. ||), with c 0.6%, which is typical of what one 
might expect given the 99% purity specification of the 
manufacturer for the sample that was used. 
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Surfactants and water have been introduced to the 
sample in order to make the emulsion droplets, and there- 
fore it is reasonable to expect additional impurities to 
end up in the emulsion droplets. A similar fit to the 
melting curve of emulsion droplets shown in Fig. ^ shows 
an increase in the impurity concentration to about 3.3%. 
Given the high purity and small solubility of water in 
C16, the surfactants are likely to be the main source of 
these additional impurities. 

Given that there is a range of melting temperatures, 
which we attribute to varying impurity concentrations in 
the droplets, we propose the following scenario in order 
to explain the observed behavior: The system is presum- 
ably in thermodynamic equilibrium prior to the initia- 
tion of nucleation events. It is reasonable to suppose 
that when a droplet freezes, some of the impurities in- 
side or at the surface of the droplet are expelled to the 
surrounding water, since the impurities are less soluble in 
solid than in liquid C16. These impurities consequently 
redistribute themselves among the water and remaining 
liquid C16 droplets in order to restore thermodynamic 
equilibrium. This happens by diffusion through the con- 
tinuous phase, which occurs rapidly on the time scale 
over which the fraction of frozen droplets varies. Some 
fraction of these impurities, determined by their relative 
solubility in water vs. C16, end up in the remaining 
liquid droplets, whose enhanced impurity concentration 
decreases their effective undercooling and reduces their 
nucleation rate. Thus, the first droplets to freeze have 
a higher melting temperature than the last droplets to 
freeze since they have a lower impurity concentration, 
which gives rise to the non-exponential behavior during 
isothermal nucleation. 

The theoretical analysis can easily be extended to ac- 
count for such a mechanism, which can then be compared 
with the experiments to test the hypothesis. We do this 
in the remainder of this section. 



B. Return to Isothermal Nucleation 

The leading behavior manifested by an impurity me- 
diated suppression of the undercooling can be expressed 
as a linear decrease in the relevant melting temperature 
as a function of solid C16 fraction n. Thus, the time 
evolution of n at a fixed temperature T is given by 



dn 
~db 



(1 - n) AV exp ■ 



9. 



(T m — T — fin) 



(19) 



where T m is the initial melting temperature of the 
droplets and j3 is the rate at which the melting point 
of subsequent droplets is reduced. In general, (3 depends 
on the volume fraction of the droplets in the emulsion 
and the relative solubilities of the impurities involved in 
the mechanism in water and C16. Being unable to de- 
termine independently, we use it as an additional fit- 
ting parameter to account for the experimental data. A 
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FIG. 6. Solid fraction as a function of time for isother- 
mal conditions. Using the droplet interaction hypothesis dis- 
cussed in the text we are able to reasonably fit (solid lines) 
the isothermal data from a range of undercooling values (open 
symbols) . 



least-squares fit of the numerical integration of Eq.(|l9|) 
with respect to (3, 7, and A (A was allowed to vary 
from the intial value as predicted using Eq.(|^) and the 
numerical constants in Table ||) was performed on the 
data. The comparison of experiment and theory for four 
different undercoolings is shown in Fig. ^. Compared 
to uniform nucleation or fixed-in-time heterogeneity as- 
sumptions, the agreement with experiment is excellent. 
Furthermore, the obtained value for [3 is physically very 
reasonable. The best-fit value for j3 was 0.55 K, i.e., the 
melting point of the last droplets to nucleate have been 
reduced by about half a degree. The resulting value of 
the pre-exponcntial factor was a factor of 10 -142 times 
the classical expression. This gave a value of A = 2 x 10 34 
m~ 3 s _1 at an undercooling of 15.1 K. The resulting value 
for the kinetic barrier was 7 = 9.8 ± 0.2 mJ/m 2 , in good 
agreement with tha t obta ined from the linear cooling ex- 
periments (see Sec. IV B| ). 

Actual values for the rates of transformation for the 

at 



isothermal data in Fig. @ range 1 
AT = 14.9 K to 7 x lfF 5 m" 3 s 



at AT = 15.3 K for 
early times. These rates slow down due to the impurity 
effect by a factor of « 10 after « 60 % conversion of the 
sample. 

It may be noted that the freezing curve with the fastest 
time scale (largest undercooling) lies somewhat below the 
theoretical prediction. We attribute this discrepancy to 
the rapid liberation of latent heat upon the freezing of the 
droplets, which raises the sample temperature and slows 
the freezing of the remaining liquid droplets. Larger un- 
dercoolings (not shown) all seem to suffer increasingly 
from this effect. In Fig. o we have also adjusted the melt- 
ing temperature of C16 from the value given in Table | of 
T m0 = 291.32 K to T m = 290.95 K. This downward shift 
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of the melting temperature from that of the pure material 
reflects the initial impurity concentration in the emulsion 
droplets prior to the freezing-induced enhancement, as 
well as the accuracy in the calibration of thermistors used 
to measure sample temperature. (Note that the melting 
temperatures for the isothermal and linear cooling ex- 
periments were measured with different thermistors for 
emulsion samples prepared at different times.) 



C. Return to Linear Cooling 

We may further test our hypothesis by revisiting the 
linear cooling problem. If the transformation rate is slow 
compared to the time over which impurity molecules are 
able to diffuse from one emulsion droplet to the next, we 
may assume that the impurities remain in equilibrium 
over the course of the experiment. We can estimate this 
equilibration time as: Tdm ~ R 2 /D [(1 — n) r/>] 2 ^ 3 where 
R is the radius of an emulsion droplet, D is the diffu- 
sion coefficient for an impurity molecule in water (which 
we may take to be on the order of a typical molecular 
diffusivity) and <p is the volume fraction of the emulsion 
droplets. We then find that the assumption of impurity 
concentration equilibration remains valid up to transfor- 
mation rates, dn/dt, on the order of r^g ~ I0 2 Hz. The 
experimental maximum transformation rate (taken from 
the maximum slope of the data in Fig. ^) is on the order 
of 10 _4 Hz <C T^g. Clearly this estimate applies only to 
the situation where the total conversion to the solid is 
incomplete (as is the case in our isothermal experiments 
- see Fig. ^) since the diffusion time to another unfrozen 
droplet becomes infinite as the fraction of remaining liq- 
uid droplets (1 — n) drops to zero. In the linear cooling 
experiments where the transformation is complete, the 
assumption of impurity equilibration holds until less than 
0.1% of the liquid droplets remain, so the observation of 
the breakdown of this assumption is experimentally in- 
accessible. 

Upon modifying Eq. (|^) to account for the suppression 
of the melting temperature, time evolution of n for linear 
cooling is given by 



dn 
~dl 



— = (l-n)AVe T^+i^y 2 



(20) 



The integration of Eq. (J2(j) can be performed numerically. 
The result of this integration along with the linear cooling 
data is shown in Fig. |?] for two different cooling rates. 
Once again, [3 is used as a fitting parameter, with a best 
fit value of 0.18 K. Values of A and 7 are as given in 



Sec. pyq . 

Comparing Fig. |^ with Fig. |] one immediately notes 
that the discrepancy in maximum slope between the un- 
modified theory and the data has now disappeared. The 
difference in the values of [3 between the linear cooling 
and isothermal quench experiments is not unexpected, 
given the fact that the value of [3 is sample dependent and 
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FIG. 7. Solid fraction as a function of temperature at two 
different linear cooling rates (open symbols) and the predic- 
tion produced by the numerical integration of Eq. ^ (solid 
lines). 



different emulsion samples were used for the two types 
of experiments. A more detailed analysis of the impu- 
rity mediated interaction hypothesis, including a better 
quantitative analysis of /3 and further implications of this 
hypothesis on other types of experiments is the subject 
of ongoing work |l|. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we report on the results of nucleation 
studies in C16 emulsions during isothermal quench and 
linear cooling conditions. By application of recently 
developed emulsion preparation techniques to produce 
nearly monodisperse emulsion samples and use of x-ray 
scattering for direct observation of the crystalline phase, 
we greatly reduced many of the uncertainties in the in- 
terpretation of experimental data. In doing so, we have 
found potentially important effects not accounted for by 
the standard interpretation of uniform nucleation. In 
accord with previous measurements of this sort by Turn- 
bull and Cormia H, the isothermal quench experiments 
do not show the standard exponential growth and satu- 
ration of the fraction of solid droplets over time. Instead 
we observe a substantial decrease in the rate of nucleation 
events at later times, when a significant fraction of the 
system has already transformed to the solid phase. In 
addition, we find that in the linear cooling experiments 
the maximum transformation rate predicted by the sim- 
ple nucleation theory is larger than that of the data. 

Such discrepancies have previously been attributed to 
the existence of some sort of temporally fixed heterogene- 
ity in the emulsion droplet population. One proposed 
source of heterogeneity has been the volume polydisper- 
sity in the emulsion droplets. However, using recently 
developed emulsification techniques, we can better con- 
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trol the droplet polydispersity and still not recover sin- 
gle exponential behavior. Furthermore, by attempting to 
fit isothermal data at different temperatures to a given 
droplet size distribution, we find that no single distribu- 
tion is capable of fitting the data at different levels of 
undercooling. Along similar lines one may attempt to 
postulate a distribution of the number of heterogeneous 
nucleators in the droplets. The best attempts to fit the 
data with such a scheme result in unaccountably broad 
distributions of impurity concentrations in the emulsion 
droplets. Such distributions are highly unlikely in light 
of the central limit theorem. 

In contradistinction to both of these ideas, we propose 
that the observed heterogeneity in the emulsion droplets 
is dynamically generated by the nucleation process it- 
self. Freezing droplets expel impurities (some of which 
might be attributed to the addition of surfactant) which 
then go into solution in the remaining liquid C16 droplets 
in order to restore thermodynamic equilibrium. By do- 
ing so, these mobile impurities reduce the melting tem- 
perature of the remaining liquid droplets, which experi- 
ence a reduced thermodynamic driving force towards so- 
lidification and consequently a reduced nucleation rate. 
This proposal is able to account for both the isother- 
mal and linear cooling data at a variety of undercoolings 
and cooling rates, respectively. To a good approxima- 
tion, the magnitude of this effect can be parametrized 
by a sample dependent parameter (3, which expresses 
the rate at which the melting point of remaining liquid 
droplets is suppressed as a function of the solid droplet 
fraction. To obtain accurate quantitive measurements 
of either the activation barrier toward nucleation or the 
pre-exponential factor, this effect needs to be taken into 
account. Note that measurements described in this pa- 
per were performed on emulsion samples prepared with 
a mixture of anionic surfactants. Measurements of a C16 
emulsion sample prepared with a single nonionic surfac- 
tant does not appear to show such a large impurity effect 

Further tests of this theory can be anticipated. It 
should be possible to directly detect and measure the 
suppression of melting temperature by an experiment 
where the emulsion is partially melted and subsequently 
refrozen at a temperature where no new nucleation is al- 
lowed. Such careful studies of the remelting process cou- 
pled to a more detailed theory for the impurity expulsion 
and transport should lead to a better understanding of 
whether and how the dynamics of the freezing process 
effects the results of emulsion nucleation studies. Addi- 
tional work along these lines is currently in progress [[H) . 
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